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The Bures metric, defined on the nondegenerate density matrices, has been the object of considerable study [jl]-[|. 
It is the minimal member of the nondenumerable family of monotone metrics HQ] . Other members of this family of 
particular note are the "Bogoliubov-Kubo-Mori" (BKM) metric the maximal monotone metric, as well as the 

"Morozova-Chentsov" JTcj ] and "quasi-Bures" pd[ | ones, the last yielding the minimax/maximin asymptotic redundancy 
in universal quantum coding [ ]l2| . Interestingly, these (operator) monotone metrics correspond in a direct fashion to 
certain "measures of central tendency" , with, for example, the Bures metric corresponding to the arithmetic mean, 
(x + y)/2, of numbers x and y BJtJ. All these monotone metrics constitute various extensions to the quantum domain 
of the (unique) Fisher information metric in the classical realm where the objects of study are probability distributions 
(rather than density matrices) Jl3| (cf. JlJ]). However, the only one of these monotone metrics that can be extended 
to the boundary of pure states yielding the standard Fubini-Study metric on this boundary is the Bures (minimal 
monotone) one || sec. IV]. 

Dittmann has shown that "the connection form (gauge field) related to the generalization of the Berry phase to the 
mixed states proposed by Uhlmann satisfies the source-free Yang-Mills equation *D * Deo, where the Hodge operator 
is taken with respect to the Bures metric on the space of finite-dimensional density matrices" |l5| (cf. |L6]]). (Let us 
also indicate here that the BKM metric has been shown to be the unique monotone Riemannian metric for which 
two "natural flat" connections — the "exponential" and "mixture" ones — are dual ||.) Here we report our efforts 
to evaluate — using numerical (lattice-gauge-like [|l7j ) methods — the action functional for this Yang-Mills field over 
the eight-dimensional convex set (M) equipped with the Bures metric of 3 x 3 density matrices, as well as for the 
two self-dual components of this field. (Using an exact version for the 2x2 density matrices of the methodology 
described below, we have found a value of precisely three for the action of the corresponding Yang-Mills field. The 
Chern-Simons functional is zero.) 

To proceed, we exploit our recent work fll8| , p"9| in determining the elements of the Bures metric using a certain 
Euler angle parameterization of the 3x3 density matrices p0| . (Earlier, Dittmann had noted jl|] "that in affinc 
coordinates (e. g. using the Pauli matrices for n — 2) the [Bures] metric becomes very complicated for n > 2 and no 
good parameterization seems to be available for general n".) This allows us, among other things, to employ as our 
parameter space four (noncontiguous) eight-dimensional hyperrectangles (each having three sides of length 7r, three 
sides of length one j in length and one cos -1 -^j), rather than less analytically convenient ones, such as that 
discussed by Bloore 0, Fig. 3], involving spheroids and "tetrapaks". (The four hyperrectangles are not adjacent due 
to the fact that two of the six Euler angles used in the parameterization have disconnected ranges. This situation — 
which results in different normalization factors only — has only come to our attention relatively recently, and serves 
as an erratum to earlier analyses |2^,|23| , which had relied upon work of Marinov |2^] , without taking into account a 
correction || to §| (cf. @).) 

To be specific, we seek to implement the formula |27], (1.1.10)] for the Yang-Mills functional YM(A) of the appro- 
priate connection A 



r W( 11 = i / \ p M 2 dv g . (1) 
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Here Fa is the curvature of A and, for our purposes, dV g is the volume form of the Bures metric g on the eight- 
dimensional manifold M, composed of the three-level quantum systems. More detailedly |2j] ], 

|F A | 2 = <F aij ,F akl >g ik 9 j \ (2) 

where (g-ij) is the metric tensor of g and (g 13 ) is its inverse. (Formulas for V g , (gij) and (g lJ ) 1 using the Euler-angle 
parameterization pOJ , have been derived in |Q (cf. p[ . Let us indicate here that since we directly adopt formulas 
from a number of different sources, our overall notation may not be as uniform as fully desirable, though hopefully 
readily enough comprehensible to the informed reader.) 

To achieve the decomposition Fa = + F^ of Fa into the sum of "anti-self-dual" (Fj[) and "self-dual" (-F7) 
parts eq. (2.9)] (29| eq. (3.4)], we use the formulas first given by Corrigan et al and presented again numerous 
times (e. g. |3^-Q). The components of the anti-self-dual field satisfy a set of seven linear equations (corresponding 
to an eigenvalue 1) ||(], eq. (3.39)] and those of the self-dual field, a set of twenty-one linear equations (corresponding 
to an eigenvalue of -3) J|(], eq. (3.40)]. The set of seven equations has been viewed as more fundamental in nature 
than the set of twenty-one equations by Figueroa-O'Farrill, who is also somewhat critical of the application of the 
terms "self-dual" and "anti-self-dual" in this octonionic setting |n| . 

For the required connection form A, Dittmann has presented the general formula Jl5|, eq. (9)] 

A= . , (W*T -T*W) (3) 
L + R 

Here the elements of T lie in the tangent space to the principal f/(7i)-bundle, the manifold of invertible normalized 
(TrW^*W = 1) Hilbert-Schmidt operators, while L is the operator (depending on W) of left multiplication by the 
density matrix p = WW* and R the corresponding operator of right multiplication. Also, L and R are the counterpart 
operators for p = W*W . Since the Euler-angle parameterization of the 3x3 density matrices presented in |^0| is of 
the ("Schur-Schatten") form, 

p = UDU*, (4) 

U being unitary, U* its conjugate transpose, and D the diagonal matrix composed of the three eigenvalues of p, one 
can immediately express W as UD 1 / 2 !/*, with W = W*. 

In the non-trivial task of implementing formula (||) for the three-level quantum systems, we relied upon the implicit 
relation (equivalence) between two formulas for the Bures metric for n-level quantum systems |]35| , eqs. (2) and (16)], 

9 = l TTdp LTR dp > (5) 

and 

1 " 

g = - J £a ij Trdpp i - l dppi-\ (6) 

(The somewhat involved formula for the coefficients a^, functions of elementary invariants, is stated in Proposition 3 
of |35|| .) In other words, we used the formula for the connection 

1 3 

A = -J2a ij p i - 1 S(P- 1 ) (7) 

where S — W*T — T*W. (We were able to express some of the individual components of A in relatively simple forms, 
but others remain, at this stage, quite cumbersome in nature.) 

Let us note that in a previous similarly-motivated study ]l9| to that here, we pursued many of the same questions, 
but somewhat misguidedly employed for our analyses the curvature (R) of the Bures Riemannian metric, rather than 
the curvature (Fa) of the gauge field A, so the interpretation to be given to the results there, in retrospect, is not 
altogether clear. (By embedding the spin connection in the gauge connection, one can construct a self-dual gauge 
field directly from a self-dual metric pq].) In that analysis M], numerical evidence led us to conclude that for the 
8x8 skew-symmetric matrices of the R a bcd, holding c, d constant, one of the four pairs of imaginary eigenvalues is 
always degenerate, that is (0,0) cf. p^gq|. 
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In p9[, we pursued both a Monte Carlo numerical integration scheme, randomly selecting points lying in the ranges 
of the eight parameters, as well as a scheme in which the points were systematically selected as the grid points of a 
hyperrectangular lattice. Here we have chosen to focus on a lattice scheme, employing 1,024 = 4 x 256 grid points 
(cf. ]l7[|). The coordinates of the nodes of the lattice were chosen in the natural (symmetric) manner, so as to divide 
each individual side of the four hyperrectangles into three intervals, the middle interval being twice the length of the 
two end ones, which under wrapping (identifying the end points) form a single interval equal in length to the middle 
one. Based upon this lattice scheme, we approximated the value of the Yang-Mills functional (|l|) for the Yang-Mills 
field F A itself to be 137.653, for the "anti-self-dual" part F| to be 222.806, and for the "self-dual" part F7 to be 
176.194 (Table I). 



quantity 


F A 


FX 


FX 


Yang-Mills functional 


137.653 


222.806 


176.194 




37.2769 


15.0131 


1.1757 




162.468 


97.6351 


4.77715 




129.285 


72.9629 


8.54855 


H A 


-1.35 ■ 1CT 11 


4.15623 


35.6896 


H 2 


-1.74 ■ KT 11 


5.4683 


71.7151 


Hi - Hi jl 


-4.77 ■ lO -1 ^ 


1.42208 


-.167908 



TABLE I. Lattice approximations, based on 1,024 = 4 x 2 8 grid points, to various quantities of interest for the 
eight-dimensional Yang-Mills field Fa over the three-level quantum systems, and for its anti-self-dual (F£) and self-dual 
(F^) constituents 
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The explanation for the Yang-Mills functionals being larger for both F A and F A than for Fa itself, in apparent 
contradiction to the general assertion in |2^, sec. 2], may be that the topological types of F A and F A are different 
from that of Fa- (I thank Y.-H. Gao for suggesting this possibility.) It should be noted, though, that in more than 
four dimensions, the Yang-Mills functional itself need not be a minimum of the action sec. 2.4] pp| . 

Motivated by the work of Bilge, Dereli and Kogak (BDK) |34|] concerning >S'0(n)-bundles in eight-dimensional 
spaces, we computed the integrals over M of (Fa, Fa) 2 , (F\, F\) and (tiF A , trF|), and also similarly for the self-dual 
and anti-self-dual parts of Fa ■ BDK considered a generic action density of the form 

a(F A ,F A ) 2 + b(Fl,Fl)+c(trFl,tvFl). (8) 

(The action density proposed by Grossman, Kephardt and Stasheff [[ll| corresponds to the choice a = c = 0.) BDK 
derived a topological bound 

/ (F A ,F A ) 2 >k f Pl (E) 2 (9) 

J M JM 

that is achieved by "strongly (anti)self-dual" Yang-Mills fields (those conforming to the set of twenty-one — not the 
set of seven — equations of Corrigan et al), where pi is the first Pontrjagin class of the SO(n) Yang-Mills bundle E 
and k is a constant. This is to be compared with a topological bound obtained earlier [^dj 

/ (F 2 A,F 2 A )>k' f P2 (E). (10) 

JM JM 



We also approximated for the three fields (Fa,F\ and F A ) the quantities |l2], eq. (Al)], where (using F to 
' A andF+; 



generically represent Fa, F a and F A ) 



and 



/ ^^Lypar\Kri^(F^uFp a F T \F Kir ^ (1 1) 

JM 

^tivparXK'q^(Fp J jjFp a ) tr (F T \F Kr ^ ) , (1^) 

JM 

for which H4 — H 2 /2 is proportional to a certain topological invariant, qs (4^, eq. (Al)] (e^ pcrT \ Kri is the completely 
antisymmetric Levi-Civita symbol). Let us note that Table I indicates that for the Yang-Mills field Fa, the terms 
i?4 and H 2 are both essentially zero. (Hi is proportional to the i-th Chern number. "The Chern numbers are those 
obtained by integrating characteristic polynomials of degree dim M over the entire manifold M" ^3|, Def. (2.5.2)].) 

In contrast to the Yang-Mills functional ([[]), we see from Table I that (following the work of BDK [f34| ) the 
integrals over the Ricmannian manifold composed of the three-level quantum systems endowed with the Bures metric 
of (F, F) 2 , (F 2 ,F 2 ) and (tri^ 2 ,trF 2 ) are largest for the Yang-Mills field F A H), smaller for the anti-self-dual part 
F A and smallest for the self-dual component F A . 

In the context of 5 l O(n)-bundles and 2-forms in 2m dimensions represented by 2m x 2m skew-symmetric real 
matrices, BDK derived the result that 

(F A ,F A ) 2 >l(trF 2 ,tvF 2 ). (13) 

We see from Table I that for our SU (3)-bundle, on the other hand, this result does not appear to hold for all three 
fields. 

For analyses of non-self-dual Yang-Mills fields (in four dimensions) see p5|-^8| and (more generally, in 2n dimensions) 
© K 11). 

The computational demands in deriving (with the use of MATHEMATICA) the results reported above in Table I 
have been considerable — due, in particular, to the still somewhat cumbersome formulas for certain of the elements of 
the Bures metric and the term W*T — T*W [occurring in the formula (||) for the connection A] that we have derived 
Hpl, using the parameterization in Q. To employ instead of the 1, 024 = 4 x 2 8 point lattice, a 26, 244 = 4 x 3 s 
one, for example, presently seems too much for us to reasonably accomplish. (Actually, it appears that since the 
results for the four noncontiguous hyperrectangles are all identical, we would only have to undertake 6, 561 = 3 s 
independent evaluations.) However, we have been able to refine the lattice in two of the eight directions, so that we 
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deal with 2, 304 = 4 x 3 2 x 2 6 points. If we choose these two directions (more or less arbitrarily) to correspond to one 
of the six Euler angles (a in the notation of Jl8],[l9]]) and one of the spherical angles parameterizing the eigenvalues 
(#i), then we obtain the outcomes reported in Table II. 



quantity 


Fa 


n 


FX 


Yang-Mills functional 


197.458 


324.53 


256.239 




37.9854 


15.3356 


1.20481 




162.107 


98.5983 


4.92389 


f„(trF*,trF*) 


132.557 


74.8023 


8.77179 


Hi 


-3.08- 10" 12 


4.20952 


36.6592 


H 2 


4.13 • 10 _ii 


5.28634 


73.5915 


H 4 - H 2 /2 


-2.37 • 10 _ii 


1.56635 


-.136522 



TABLE II. Lattice approximations, based on 2, 304 = 4 x 3 2 x 2 6 grid points, to various quantities of interest for the 



eight-dimensional Yang-Mills field Fa over the three- level quantum systems, and for its anti-self-dual (Fj[) and self-dual (F A ) 
constitutents 
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We see (quantitative, not really qualitative) changes in Table II from the values reported in Table I, witt quite 
substantial increases in the case of the Yang-Mills functionals, and comparatively modest differences in the other sets 
of indices Q . 
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